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Abstract. Cooling of the mechanical motion of a GaAs nano-membrane using the 
photothermal effect mediated by excitons was recently demonstrated by some of us 
[K. Usami, et al, Nature Phys. 8, 168 (2012)] and provides a clear example of the use 
of thermal forces to cool down mechanical motion. Here, we report on a single-free- 
parameter theoretical model to explain the results of this experiment which matches 
the experimental data remarkably well. 
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Figure 1. Schematic of the experiment. A Fabry-Perot cavity is bounded by a perfect 
mirror on one end and a GaAs membrane on the other. Drumhead vibrations of the 
membrane, indicated by the dashed lines, couple to the cavity field and to the excitons 
inside the membrane. In the experiment, L = 2.9 cm, d = 160 nm, and the membrane 
is not perfectly flat; this latter feature is important to the mechanism we consider, as 
explained in Sec. 4. 

1. Introduction and motivation 

Micro- and nanomechanical systems in the quantum regime offer exciting perspectives 
for fundamental tests of quantum physics as well as for quantum technological 
applications. It is a unique feature of these systems that they can be strongly coupled 
to a plethora of other quantum systems: optomechanics explores interactions with light 
based on radiation pressure [1, 2], dipole gradient [3, 4], or photothermal forces [5, 6, 7]; 
electromechanics investigates the coupling of electronic and mechanical degrees of 
freedom [8, 9]; and magnetic forces can couple a mechanical oscillator to magnetic 
moments [10], even of single electrons [11]. It is this versatile nature of mechanical 
systems which makes them attractive as basic building blocks for hybrid quantum 
systems [12]. In the current paper we elaborate on the recent findings reported in 
Ref. [13] demonstrating the interplay of photonic and electronic degrees of freedom in a 
micromechanical semiconductor membrane. 

The photothermal effect [5, 6, 7, 14, 15] elegantly overcomes one fundamental limit 
encountered by any optomechanical cooling mechanism based on the radiation-pressure 
interaction [1, 2]. That is to say, a single photon of frequency Co>l can only provide an 
energy change AE ~ (v/c)hu^ upon reflection off a mirror with velocity v, and this due 
to the Doppler shift [7], but absorption of the same photon, as per the photothermal 
effect, implies AE ~ hu^. How this AE translates to a change in motional energy is a 
less well-defined concept, but it is clear that, in principle, the latter effect can give rise 
to cooling forces that eclipse the radiation-pressure force produced by the same number 
of photons. 

In the semiconductor GaAs, the absorption of an above-bandgap optical photon and 
subsequent decay of the associated bound states [16, 17] is a complicated process; the 
energy liberated by the photon first creates an electron-hole bound pair (an exciton), 
which decays by scattering phonons throughout the structure of the material. This 



Exciton-mediated photothermal cooling in GaAs membranes 



3 



(a) 




(b) 




radiation pressure 



\ exciton creation 
fi CJ f2 in \ and radiative decay 



deformation potential 



Figure 2. (a) The full set of interactions that our model describes; each interaction 
is labelled by the relevant coupling constant in the effective description of Sec. 3. 
(b) Adiabatic elimination allows us to describe the system in terms of one cycle of 
interactions, with the mechanics modifying the optical fields, which act on the excitonic 
fields, which in turn couple to the mechanics. A memory kernel, used to describe a 
delayed interaction, will later be introduced into the segment connecting c\ and b. 



scattering process, which manifests itself primarily as the transport of heat and takes 
place over the thermalisation time r t h, changes the properties of the material. The 
drumhead modes of a membrane are critically dependent on these properties, and this 
process therefore couples the absorption of the light to the motion of these modes. 
Our main aim in this paper is to describe a phenomenological Hamiltonian model for 
this process and to predict the cooling or heating effect imparted by this photothermal 
interaction. We take particular care to model the details of the experimental system of 
Ref. [13]; we notice, for example, that the coupling of the excitonic field to the continuum 
of input field modes is crucial to describing the experimental data. 
This paper is structured as follows. In the next section we shall build our system 
Hamiltonian and derive the equations of motion of the field operators. Following this, 
we shall insert a memory kernel in these equations to account for the delay in the 
thermalisation process. After we fit our model to the experimental data, we briefly 
discuss the possibility of using pure deformation-potential effects to achieve cooling in 
similar systems, and then conclude. 

2. System Hamiltonian and Equations of motion 

Let us start by considering the physical system represented in Fig. 1: a Fabry-Perot 
cavity is bounded by a perfect mirror on one end and a semi-transmissive GaAs 
membrane on the other. We shall describe the exciton fields inside the membrane 
by means of a bosonic approximation, which is valid when the exciton population is not 
too large, and assign to these fields the operators c\ (frequency uj\ and decoherence rate 
7a), where A is some, possible continuous, index labelling the modes. Any sum over A 
is to be interpreted as either a sum, if A is discrete, or an integral, otherwise. 

These fields interact with the cavity field, a (frequency u c and HWHM linewidth 
k c ), as well as with an infinity of field modes that represent the free field that 
forms the input to the cavity field. The excitonic fields also interact with our chosen 
mechanical mode, which we describe using the operator b, and which has a mechanical 
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frequency u m and HWHM linewidth K m . 

The Hamiltonian H is made up of four different contributions. The free Hamiltonian 
reads (we shall take h = 1 throughout this paper for conciseness of notation) 

-fffree = u c a j a + duj wa+ a w + ^ w a4<?a + ujfib , (1) 

J A 

and describes the fields in the absence of any interaction or dissipation. The next set of 
terms describes dissipation, and reads 



#diss — ^ U _ d^a)) + -f/mcch.diss + -^exc,diss ] (2) 




we have chosen to write down explicitly only the Hamiltonian describing the cavity 
field decay. i^ me ch,diss (-^cxc,diss) similarly describes the dissipation of the mechanical 
(excitonic) operator (s). There are two sets of interaction terms, illustrated pictorially 
in Fig. 2(a); the first couples the excitons to the electric field: 

H ahs = duJ ^P^H" + fi c,A^4 + H - c - > ( 3 ) 

where fi c ,A ; which is assumed to be real, is the coupling rate of exciton mode A with 
the cavity field, and £l- m> \ with the free-field modes. The coupling of the exciton modes 
to the free field cannot be neglected in this case: the membrane is thick enough for 
interference effects between the input field and the cavity field to be significant, and 
|^in,A/^c,A| to be of order unity. The two terms in H^s therefore interfere, leading to an 
asymmetric cooling spectrum (see Sec. 4, below, and Fig. 3); this phenomenon is closely 
linked to the Fano line-shapes observed in optomechanical systems where the mechanical 
oscillator is coupled both to a cavity field and to the free field [18, 19], although the 
dominant coupling of the optical fields to the mechanics is indirect in the present case. 
The second set of interaction terms describes the coupling of the mechanical motion to 
the cavity field and to the excitons [20, 21]: 

# mcch = g tfa(b + S f ) + ^a,a4ca' (b + P). (4) 

A, A' 

go is the usual (radiation-pressure) coupling constant and h\ t \i describes the 
deformation-potential coupling. We omit any direct coupling of the motion to the free 
field, since such effects would be very small [19] compared to the terms in the preceding 
equation. Finally, we can write 

H = H iree + H diss + if abs + H mech . (5) 

This Hamiltonian can be used to generate the Heisenberg-Langevin equations of motion 
for the field operators. We choose to work in a frame rotating at the frequency of the 
driving field, cjl, and define the detunings A c = wl — w c and A a = wl — u)\. At this 
point we shall make two further assumptions regarding the exciton fields. The quantities 
Q C) x Q c , r2 in A —7- Q in , h\ \i — > ho, and 7a — > 7 are assumed to be independent of the 
index A. We also assume that the exciton density of states is constant in the relevant 
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region, which is a good approximation for pumping well above the band-gap energy. In 
order to avoid introducing new symbols we shall displace each operator 6 by its mean 
value o: 6 — > o + 6, where the operator on the right-hand side has zero mean. Therefore, 
all operators in the following will have zero mean. The linearised equations of motion 
can now be written as 

a = — (k c — iA c )a — y/2K^d- m — igo<x{b + ¥) , (6) 

b = — (ft m + iuj m )b — \f2K m i> m — ig (a*a + da 1 ) — i ^(/i*c A + hc\) (7) 

A 

and 

ca = -(7 - ^A A )c A - v^cx.in - ih{b + b ] ) - i(Q in + Q c )d 

^in^in j (8) 

where terms without a significant contribution were dropped, and where we defined 
h = h J2\C\-- The operators c Ajin describe the zero-mean Langevin forces associated 
with the excitonic modes, whereas the optical field input operator a in is defined as per 
the usual input-output theory (cf. Ref. [22, §5.3]), but we note that the interaction 
of the excitons with the free field modifies the input-output relation for the system, 
yielding 



y c A 



«out = a in + V2fi^a - i\ — fi* n \ c A . (9) 



The mean values of the fields satisfy (c in = 0) 

— (k c — iA c )d — i(Q c — Qi n )* c A — A/2K c a in = , and (10) 

A 

-(7 - «A A )C A - i(fl in + ttc)d - ij — fiinOin = . (11) 

V «c 

We have absorbed b into an effective redefinition of A c and A A , and therefore we have 
6 = 0. It is now apparent that the equations for the optical and mechanical fields involve 
only sums of the type c\ or c~\. These sums can be performed easily due to our 
assumptions, yielding 

2__ \ K c 7 + vQ in (Q c - ^in)* / 12 ) 

and 



ca- ^ u []j Kc a ^J ^ c _ iAc j J + u ^ c + n . n ^ c _ Q . u y^ y ) 

where v = 7^ A (7 — iA A ) _1 6 IR + accounts for the number of exciton modes we are 
interacting with. 



Exciton-mediated photothermal cooling in GaAs membranes 



6 




-1 -0.5 0.5 1 1.5 2 -1 -0.5 0.5 1 1.5 2 

Cavity detuning, A c / (2re c ) Cavity detuning, A c / (2k c ) 

(c) Dataset 3; Al = 852 nm (d) Dataset 4; Al = 852 nm 



Figure 3. Fits to four sets of experimental data; Dataset 1 is the set reported in 
Ref. [13]. The four sets of data differed in the location of the driving beam on the 
membrane and in the driving wavelength Al used. All other system parameters are 
reported in the text and in Table 1. 

3. Adiabatic elimination and memory kernel 

The system we aim to describe has a hierarchy of dynamics determined by 7 ^> k c ^> 
K m ,u m . Sequential adiabatic elimination of the fields is therefore possible, first solving 
the equation of motion for the excitonic fields, on whose time-scale a is approximately 
constant, and then for the cavity field. We concentrate on the regime where the 
exciton-mediated effect dominates over the radiation pressure force. At this level of 
approximation we can describe the web of interactions in our system via the following 
process: 

(i) the position of the mechanical oscillator changes the photon number inside the 
cavity {b affects a), 

(ii) the exciton populations follow the changing cavity field (a affects d\), and finally 

(iii) the changing exciton populations modulate the mechanical properties, coupling to 
b (ca affects b). 

Fig. 2 illustrates the different interactions that take place in the system we are describing, 
together with the reduced system that results after adiabatically eliminating the optical 
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and excitonic fields. 

With the above in mind, and ignoring the contributions from the input field, we can write 
the adiabatic solution of the cavity field operator familiar from the radiation-pressure 
cooling literature [23, 24]: 

"1? 1 



a Pa —ig^a 



_k c - i(A c + uj m ) K c -i(A c -u 
which is substituted into the equation of motion for c\ to give 



g a(n c + O i: 



7 - « c - i(A x - A c 



1 -6+ ttI $ 



K c -i(A c + u m ) K c -i(A c -u n 



(14) 



.(15) 



We see from this equation that the function of the excitons at this level of approximation 
is to act as a channel for the absorbed optical energy to interact with the mechanics. 
Noting that 7 is much larger than all the other frequencies, we obtain 

v g a^(Q c + Q in 



A 



7 



k c - i(A c + u m ) 



+ — 7x i bi 



(16) 



Substitution of this solution into the equation of motion for b gives us the adiabatic 
dynamics when the excitons couple to the mechanics through the deformation potential. 
Because we want to describe the time-delayed effect of the excitons on the mechanics, 
however, we must introduce a memory kernel Ai(t) into the equation of motion 
for b. This follows the ideas outlined in related treatments of the photothermal 
effect [15, 5, 6, 7], and in our notation corresponds to setting: 

/oo 
M(t-T)J2cx(r)dT, (17) 

a 00 A 

and similarly for the term involving ho- This process changes the physical meaning 
of these terms. We emphasise that 77 no longer describes the deformation-potential 
coupling, but is a phenomenological coupling constant that describes the strength of the 
delayed interaction linking the exciton fields with the mechanical motion. Physically, 
Eq. (17) tells us that the modulation to the mechanical properties takes into account the 
entire history of the exciton fields. We choose to use an exponentially-decaying memory 
kernel [7]: 

M(t-r) = — e - {t - T)/T ^Q{t - r) , (18) 

Tth 

where 0(t) is the Heaviside step function and accounts for the causal nature of the 
memory kernel; this M.{t) leads to the same expressions as the il h(t)" chosen by Metzger 
and co-workers, cf. Ref. [15], upon integration by parts of the relevant terms. The use 
of a memory kernel in the equation of motion for b can be motivated by making use of 
an extended model that includes a bath of phonon modes which act as the intermediary 
between the excitons and b. Elimination of the these modes in db/dt naturally gives 
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Table 1. Experimental parameters and coupling constant extracted from the data. 
The datasets are numbered as per Fig. 3. 
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55 
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4.6 


3 


852 


25 


2.2 


55 


35.4 


7.6 


4 


852 


25 


2.2 


55 


35.4 


6.2 



rise to a time-integral of a sum of decaying exponentials, which we identify with an 
exponentially-decaying memory kernel; this process is outlined in Appendix A. 

A brief note about the effect of noise terms is due. It lies outside the scope of 
this paper to consider the effects of noise on the limits of this cooling mechanism; being 
interested in cooling rates in this article, we accordingly discard such terms. As discussed 
in Ref . [7] , the nature of the photothermal effect does not preclude reaching the ground 
state, even in the bad-cavity limit, despite the absorption of light in the mechanical 
oscillator. The thermal noise induced by the absorbed light can be effectively modelled 
as a Langevin force term as in, e.g., Eq. (5) of Ref. [5] or Eq. (15) of Ref. [7]. This 
Langevin force can be viewed as having its physical origins in the thermal fluctuations 
of the phononic bath described in Appendix A, which couples the exciton modes to b. 
A unique feature of our setup lies in the interference between the cavity and 
input fields, discussed after Eq. (3) above, which could, in analogy with dissipative 
optomechanics [19], lead to a situation where the effect of the noise originating from 
the optical fields cancels out, and therefore to a more efficient cooling mechanism and 
a lower base temperature. 

4. Optomechanical cooling rate 

Proceeding from the previous section along the same lines as standard optomechanical 
theory, we can now derive a simple expression for the optomechanical cooling rate due to 
this photothermal effect. Indeed, we can show that the mechanical decay rate changes 
from K m to n m + Kth., where for T tn ^> l/u m and Q c G M. 

= P in 2 9o gth g| Re f (1 + gin/go) ( A c ntjn c + IK C ) 

th hw L (k% + A 2 c )u m T th 77 \ k c - i(w m + A c ) 

(i + ni/n c )(A c n in /n c -iK c ) 

k c - i(u m - A c ) 

with P in = hw^ I cii n 1 2 being the input power coupled into the cavity (in watts), the sum 
over A was absorbed into the phenomenological coupling constant rjth = urj, and we 
have absorbed y/v into each of Q c and Q- m - m accordance with the approximations 
made during adiabatic elimination, Eq. (19) excludes higher-order terms in Q c and Qi n . 
It is worth noting that the sign and magnitude of i] th depend on the shape of the 
membrane. The membrane in the experiment has a slight curvature, such that any 
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thermal expansion has a well-defined effect on its effective position x = (b + w)/y/2\ if 
the membrane were to be flipped over, the sign of r^h would change. For a perfectly 
flat membrane, the membrane 'would not know' which way to buckle under thermal 
expansion; ??th would then be zero and other terms would be expected to dominate. 

The expression for n t ^ depends critically on fl in , whose relationship to Q c is fixed 
by the geometry of the cavity and membrane. In the good-cavity limit, which is valid 
whenever the finesse of the cavity is 3> 1, and taking into account the large refractive 
index of GaAs, we obtain 

^ = - e i(*L<*/2-2LA c / C ) sin (A; L rf/2) , (20) 

where d is the thickness of the membrane, L the length of the cavity, and &l = cjl/c fi c 
itself can be fixed by observing the fraction / a t, s of power absorbed by the membrane, 
since it can be shown that on cavity resonance (A c = 0) and for <C k c 7, 

Us = ^ • (21) 

It is worth noting the physical significance of Eq. (21): in our model Q c and Q{ n effectively 
give rise to the imaginary part of the refractive index of the membrane, conventionally 
labelled k. / a b s therefore exhibits a wavelength-dependence, as does k, cf. Table 1. 
Independent experimental measurements give us values for: k c , f^/7, u m , and r t h- go 
is fixed by the geometry and by the reflectivity of the membrane, whereas P> m , u^, and 
A c are determined by the experiment. The only independent fit parameter in Eq. (19) 
is therefore the photothermal coupling strength r/th/7- 



5. Fit to experimental data 

We shall now use K th to model four sets of experimental data. The experimental runs 
differ in the transverse location of the membrane vis-a-vis the cavity field, and the values 
of 77th obtained for the four sets are consistent with the coupling of the excitons to the 
(2, 1) drumhead mode of the membrane. 

The independently-determined parameters used in the model and common to every 
dataset were as follows: L = 2.9cm, d = 160 nm, go = 2tcx (— 5.1) Hz, n c = 27rx258MHz, 
u m = 2tt x 23.4 kHz, and r t h = 6.6 ms. Other parameters differed between datasets and 
are listed in Table 1. This table also lists the photothermal coupling rates resulting in 
the fits shown in Fig. 3. These coupling rates, when plotted as a function of the position 
of the driving beam on the membrane, are consistent with a membrane displacement 
profile matching the (2, 1) drumhead mode, cf. Fig. 4, and therefore conform to our 
expectations. In other words, a single fit parameter, corresponding to a maximal 
coupling strength ?7{£ ax /7 = 0.099, suffices to fit all the datasets if we also take into 
account the location of the driving beam on the membrane. 
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Left edge Center Right edge 

Position of driving beam along membrane 

Figure 4. (a) The coupling parameters that fit the four data sets are consistent 
with the excitation of the (2, 1) mode of the membrane [shown in (b)]; the solid curve 
represents the sinusoidal profile of the squared amplitude of this mode, (b) Vibrational 
profile of the (2, 1) mode of the membrane. 

6. Cooling through electronic stress 

Let us now consider a different system where the mechanism that provides the cooling 
force is no longer the thermal stress set up by decaying excitons but pure electronic 
stress. For a system such as the one we described above to enter this regime, the 
photothermal effect must be switched off. In the case of GaAs, this can be achieved by 
operating at a temperature of ca. 12 K [13]. Under these conditions, the excitons interact 
with the mechanics directly, and only through deformation potential, as expressed in 
the model of Sec. 2. In an equivalent picture, we may say that the effective memory 
time is zero, and the memory kernel reduces to a delta-function: Ai(t) — > S(t). 
By placing the membrane inside a cavity [25] (fi m = 0), the finesse of the cavity can be 
made significantly larger if level of absorption in the membrane is lowered, leading to a 
correspondingly smaller Q c . This can be done in GaAs by running the experiment at 
longer wavelengths, e.g., at 884 nm, where the absorption is significantly lower than at 
870 nm. 

7. Conclusion 

We have explored optomechanical cooling through the photothermal effect in a 
semiconductor membrane. Our model uses a coupling similar to the deformation- 
potential coupling but makes use of a memory kernel to model the long thermalisation 
time typical of such structures. The introduction of the memory kernel was based 
on entirely phenomenological grounds, following Refs. [15, 5, 6, 7], but we justify the 
use of an exponentially-decaying kernel by introducing additional phononic degrees of 
freedom that are then eliminated. The resulting model only has one free parameter, 
with all others being determined independently or by the geometry of the situation, and 
provides a remarkably good fit to the experimental data. 

By using different forms of the memory kernel we can also compare the different 
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physical mechanisms in promoting optomechanical cooling. Thus, for example, an 
instantaneous memory kernel Ai(t) = 5(t) reduces our description to one taking into 
account pure deformation-potential effects, similarly to what was originally envisioned 
in Ref. [26]. Such a mechanism could be important under conditions where the 
photothermal effect is cancelled out, e.g., at temperatures where the membrane 
undergoes no photothermal deformation [13]. 
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Appendix A. Indirect exciton— motion coupling, and emergence of memory 
kernel 

The model presented in Sec. 2 serves to explain the physical processes occurring in a more 
transparent manner. As far as the interaction of the c\ with h is concerned, however, 
this model does not capture the fact that the process takes place indirectly. In this 
Appendix we will examine a more detailed Hamiltonian that leads to the same effective 
equation of motion for h. We shall introduce a phononic bath of modes (oscillation 
frequency and amplitude decay rate « At ) that serve as intermediaries between c\ and 
b. Physically, the d^ account for the nonzero temperature of the lattice making up the 
membrane. The interaction terms between the three systems can be written as [20, 21] 

fc A,A>4cy(4 + d\) + l »A$// h + dJl#) , (A.l) 

X,\',fi fi,fi' 

with the &a,av and l^i being coupling frequencies whose values we shall not specify or 
calculate. Hermiticity requires that ky t \^ = k\ x , and i M / )(U = I*,. Any sums over A 
or \i may be either discrete or continuous, as the case requires. Eq. (A.l) mediates the 
interaction between the excitons and b, and therefore replaces the second term in if mec h. 
We shall now proceed to eliminate the phononic bath modes. The equation of motion 
for the dp reads 

dfi = -(/fy + iuj^dfj, - i ^2 k\,\i )fl c{cx - % y^^djb + l^^d^P) , (A. 2) 

A, A' fi> 

where we have not written down the input noise terms, expressed in terms of the 
anti-Hermitian operators d^ in the quantum Brownian- motion damping model [27], 
since these terms have no effect on the cooling rate but help to determine the lowest 
mechanical occupation number that can be achieved through this cooling mechanism. 
The last set of terms in Eq. (A. 2) leads, both directly and through the C\, to a 
renormalisation of co m and K m due to the absorbed optical power and finite temperature 
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of the lattice, and we may therefore safely ignore it, linearise the equation of motion, 
and finally write 

dp = + iuj^df, - i J^(£:a,a',^c A / + k X ',x^c x c\,) . (A.3) 

A, A' 

Formally, then, the solution for d^ is given by 

4(0 = -i f e -(«^)(*-r) Y}kx',x,AMr) + k x ,y^c y c{(r)]dr , (A.4) 

J -°° A, A' 

noting once more that we are ignoring input noise fields. The Hamiltonian above 
therefore gives the following contribution to the linearised equation of motion for b: 

b= -i 5^(W'<£'4 + V.A'djj 
= ~ l f E e " KM(t_T) {- ie "^ ( *" r) W4[^,A,^A'£A(r) + W*4tolK 

J — oo 

- * f E e " KM( *" r) { ie ^ ( *" r) V,^M'[^,A',^A'£A(r) + fcj, A , M Cv6 t A (r)]}dr 

J —oo 

= ~ l f E e ~ KM( *" r) He~^ (t - r) W'^ + c -C-]^,A lAt c A ,c A (r)dr 

»/ — oo 

-i f E e " KM(t " T) [- ie "^ (t " T) W<'+ c - c -] fc A,A' lA4 cvci(r)dr, (A.5) 

J — oo 

with the sums running over A, A', /i, and //; 'c.c' denotes the complex conjugate of the 
preceding term. We now make the formal replacement 

E e-^He-^I^ + c.c.]^,^, ^(E^a) -M(f - r) , (A.6) 

A',/i,/i' A 

where .M(£) is a causal memory kernel that we choose to have a decaying exponential 
form, and where &A',A,/i is assumed to be independent of A and A'. Finally, then, the 
contribution to the equation of motion for b is 



A J ~°° A 

-^(E^a) M(t-r)J2c[(r)dr. (A.7) 

A J -°° A 

The memory kernel in the equation of motion for b therefore arises naturally from this 
more complete, albeit still phenomenological, model. 



References 



[1] Kippenberg T J and Vahala K J 2008 Science 321 1172-1176 

[2] Aspelmeyer M, Groblacher S, Hammerer K and Kicsel N 2010 J. Opt. Soc. Am. B 27 A189-A197 
[3] Li M, Pcrnice W H P and Tang H X 2009 Phys. Rev. Lett. 103 223901 

[4] Xucreb A 2012 Optical Cooling Using the Dipole Force Springer Theses (Heidelberg: Springer) 
ISBN 9783642297144 



Exciton-mediated photothermal cooling in GaAs membranes 



13 



[5] Pinard M and Dantan A 2008 New J. Phys. 10 095012 

[6] De Liberato S, Lambert N and Nori F 2011 Phys. Rev. A 83 033809 

[7] Restrepo J, Gabelli J, Ciuti C and Favero I 2011 C.R. Phys. 12 860 

[8] O'Connell ADet al. 2010 Nature 464 697-703 

[9] Taylor J M, S0rensen A S, Marcus C M and Polzik E S 2011 Phys. Rev. Lett. 107(27) 273601 

[10] Forstner S et al. 2012 Phys. Rev. Lett. 108(12) 120801 

[11] Rugar D, Budakian R, Mamin H J and Chui B W 2004 Nature 430 329-332 

[12] Wallquist M, Hammerer K, Rabl P, Lukin M and Zollcr P 2009 Phys. Scr. 2009 014001 

[13] Usami K et al. 2012 Nat. Phys. 8 168-172 

[14] Metzger C H and Karrai K 2004 Nature 432 1002-1005 

[15] Metzger C, Favero I, Ortlieb A and Karrai K 2008 Phys. Rev. B 78 035309 

[16] Scgall B and Marian G D 1968 Phys. Rev. 171(3) 935-948 

[17] Gartner W W 1961 Phys. Rev. 122(2) 419-424 

[18] Elstc F, Girvin S M and Clerk A A 2009 Phys. Rev. Lett. 102 207209 

[19] Xuereb A, Schnabel R and Hammerer K 2011 Phys. Rev. Lett. 107(21) 213604 

[20] Toyozawa Y 1958 Prog. Theor. Phys. 20 53-81 

[21] Steyn-Ross M L and Gardiner C W 1983 Phys. Rev. A 27(1) 310-325 

[22] Gardiner C W and Zollcr P 2004 Quantum Noise 3rd ed (Springer) ISBN 3540223010 

[23] Wilson-Rae I, Nooshi N, Zwcrger W and Kippenberg T J 2007 Phys. Rev. Lett. 99 093901 

[24] Marquardt F, Chen J P, Clerk A A and Girvin S M 2007 Phys. Rev. Lett. 99 093902 

[25] Thompson J D et al. 2008 Nature 452 72-75 

[26] Wilson-Rae I, Zoller P and Imamoglu A 2004 Phys. Rev. Lett. 92(7) 075507 

[27] Giovannctti V and Vitali D 2001 Phys. Rev. A 63 023812 



